Tricritical point of lattice QCD with Wilson quarks at finite temperature and density 
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First principle study of QCD at finite temperature T and chemical potential [i is essential for 
understanding a wide range of phenomena from heavy-ion collisions to cosmology and neutron stars. 
However, in the presence of finite density, the critical behavior lattice gauge theory without species 
doubling, is unknown. At strong coupling, we examine the phase structure on the (p, T) plane, 
using Hamiltonian lattice QCD with Wilson fermions. A tricritical point is found, separating the 
first and second order chiral phase transitions. Such a tricritical point at finite T has not been found 
in previous work in the Hamiltonian formalism with Kogut-Susskind fermions or naive fermions. 
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I. INTRODUCTION 



One of the most challenging issues in particle physics 
is to study QCD in extreme conditions. Precise determi- 
nation of the QCD phase diagram on temperature T and 
chemical potential fj, plane will provide valuable infor- 
mation for quark-gluon plasma (QGP) and neutron star 
physics. 

For QCD with two massless quarks, several approxi- 
mation models (e.g. linear sigma model, Nambu-Jona- 
Lasino model, random matrix model, statistical boot- 
strap) [1] suggest the existence of a tricritical point on 
the (/x, T) plane separating the first order transition line 
at lower T and larger //, and the second order transition 
line at higher T and smaller /i. There has been a pro- 
posal [2] for experimental search for the tricritical point, 
via event-by-event fluctuations in heavy-ion collisions. 

Lattice gauge theory (LGT), proposed by Wilson [3] 
is the most reliable non-perturbative approach to QCD, 
based on first principles. Unfortunately, it suffers prob- 
lems like the complex action at finite /i and species dou- 
bling with naive fermions. 

Kogut-Susskind (KS)'s approach to lattice fermions 
thins the degrees of freedom, and preserves the remnant 
of chiral symmetry, but it doesn't completely solve the 
species doubling problem; It breaks the flavor symme- 
try as well. Wilson's approach to lattice fermions [3] 
has been extensively used in hadron spectrum calcula- 
tions as well as in QCD at finite temperature; It avoids 
the species doubling and preserves the flavor symmetry, 
but it explicitly breaks the chiral symmetry, one of the 
most important symmetries of the original theory; Non- 
perturbative fine-tuning of the bare fermion mass has to 
be done, in order to define the chiral limit. 



In Lagrangian formulation of SU(3) LGT at finite fi, 
complex action spoils numerical simulations with impor- 
tance sampling. The recent years have seen enormous 
efforts [4-6] on solving the complex action problem, and 
some very interesting information [7] on the phase dia- 
gram for QCD with KS fermions at large T and small 
\x has been obtained. Nevertheless, to precisely locate 
the tricritical point and critical line at large /i is still 
an extremely difficult task. One has also to resolve the 
contradiction between Monte Carlo (MC) simulations at 
intermediate coupling and strong coupling analysis for 1 
staggered flavor (corresponding to 4 flavors in the contin- 
uum): MC data [6] indicate that the chiral phase tran- 
sition at ju = and some finite Tc is of first order, 
while strong coupling analysis [8,9] favors the second or- 
der transition; The author of Ref. [9] studies the phase 
diagram on the (ji, T) plane, and discovers a tricritical 
point, while in MC simulation [6], only a line of first or- 
der phase transition is found. It might well be that they 
belong to different universality classes. 

QCD at large fi is of particular importance for neu- 
tron star or quark star physics. Hamiltonian formula- 
tion of LGT doesn't encounter the notorious "complex 
action problem" . Recently, we proposed a Hamiltonian 
approach to LGT with naive fermions at finite /i [10,11], 
and extended it to Wilson fermions [12]. The chiral phase 
transition at T — and some finite \ic was found to be 
of first order. In Refs. [13,14], the authors studied the 
phase diagram of QCD with KS fermions for 2 and 4 fla- 
vors and naive fermions for 4 flavors. A line of second 
order phase transition was found in both cases, but no 
tricritical point was found at any finite T. 

In this paper, we study the phase diagram using Hamil- 
tonian lattice QCD with Wilson fermions. At strong cou- 
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pling and in the non-perturbatively defined chiral limit, 
we find a tricritical point on the (p, T) plane. The rest of 
the paper is organized as follows. In Sec. II, we derive the 
effective Hamiltonian at finite T and /i. In Sec. Ill, we 
analyze the QCD phase diagram. The results are sum- 
marized in Sec. IV. 

II. EFFECTIVE HAMILTONIAN AT THE 
STRONG COUPLING 

A. The n — case 



According to Eq.(3) and following the procedure in 
Sec. II A, the role of the Hamiltonian at strong coupling 
is now played by 

H? ff = H eff »N. (5) 

The vacuum energy is the expectation value of H — [iN 
in its ground state and also the expectation value of 
H^ff in its ground state \ fl e ff), given by 

E n = (n\H - nN\n) = (Sl ef/ \H2 ff \Sl eff ). (6) 



We begin with QCD Hamiltonian with Wilson 
fcrmions at \i = on 1 dimensional continuum time and 
d = 3 dimensional spatial discretized lattice, 

H = M^faMx) 

X 

1 ±d 

+ rEE*W(7t-^ {x, k)^{x + k) 



x k=±l 



x j=l 



l^iv(c/ p + c/;-2) 



where 



M 



rd 

TTl H , 

a 



(1) 



(2) 



m, a, r and g are respectively the bare fermion mass, 
spatial lattice spacing, Wilson parameter, and bare cou- 
pling constant. U(x,k) is the gauge link variable at site 
x and direction k. The convention 7_fc = — 7^ is used. 
E°"(x) is the color-electric field at site x and direction j 
and Up is the product of gauge link variables around an 
elementary spatial plaquette. 

For 1/g 2 « 1, one can integrate out the gauge fields 
and derive the effective Hamiltonian H e ff, consisting of 
terms with two fermions and four fermions [15,16]. 



B. The T / and /x / case 

In the continuum, the grand canonical partition func- 
tion of QCD at finite T and /i is 



Z = Tr c-P {H -» N \ f3=(k B T)-\ 



(3) 



where ks is the Boltzmann constant and ./V is particle 
number operator 



N = ^V t (z)V>(z)- 



(4) 



C. Results for large iV c 

As shown in Ref. [12], under the mean-field approxima- 
tion, i.e., by Wick-contracting a pair of fermion fields in 
the four fermion terms in H e f f , one can obtain a bilinear 
Hamiltonian in leading order of 1 /N c 

H2 ff ~H& FA =A 1 £ii(xMx) 

X 

+ (B-n)J2^(xMx)+C, (7) 



where 



A = M 



Kd 
2aN t 



(l-r 2 )v, 



Kd(l + r 2 ) ( J \ 
a \2N C J 



C=- 



Kd 

AaN c 



((l + rV 2 -(l-rV) N s Nf. (8) 



The coefficient A plays the role of dynamical mass of 
quark, v and are respectively the expectation of 'ipip 
and t^ip in |fi e //)j divided by N s and Nf, i.e., the total 
number of lattice sites and the number of flavors. 



K = 



1 



9 2 C N 



(9) 



is the effective four fermion coupling constant. Here 
Cn = (N 2 — 1)/(2N C ) is the Casimir invariant of the 
SU(iV c ) gauge group. For N f = 1,2 and N c > 3, the 
next-to-leading corrections are small. 

In the effective Hamiltonian (7), there are three input 
parameters: r, m and \i. Suppose we study the phase 
structure of the system in the chiral limit, reached by 
fine-tuning the bare quark mass m so that the pion be- 
comes massless. The parameter M in the chiral limit is 
given by [12] 



M M chira i = - 



Kdr 2 
aN r 



-v. 



(10) 



In the chiral limit, there are only two free parameters 
left: r and fi in Eq. (7). The vacuum energy is 



2 



/ Kdr 2 
E n = 2N c NfN s I M chiral (n + n - 1) — (n - n + 1) 

Kdr 2 , n ^ Kd , 

H (n +h +1 — 2n) H (n + n - 2nn - 1) 

a 'a 



(11) 



- n(n - n + 1) , 



where the quark number n and anti-quark number n are 
constrained in the range of [0, 1] and determined by min- 
imizing the grand thermodynamic potential 



Here 



$ = E n - TS. 



S = _ 2N ° N f N » [ nlnn + (1 - n) ln(l - n) 



(12) 



+ n In n + (1 — n) ln(l — n) 



(13) 



is the entropy. 

Once n and n are known, the chiral condensate and 
quark number density in the leading order of l/N c can 
directly be obtained 

(vV) _ <"IE,g(»W»)l") : , 



JV>JV a 
= 27V c (n + n- 1), 



2N c NfN a 



2A f 



= n — n. 



(14) 



According to Eqs. (8) and (10), in the chiral limit, the 
dynamical mass of quark m^n = A is proportional to ?} 
or the chiral condensate. For convenience, we rescale the 
chemical potential and temperature as [t! = \ij{Kdja) 
andT' = T/(Kd/a). 

Figures 1, 2 and 3 show the 3D plot of $ as a func- 
tion of n and n for r = 1 and some (//,T'). At 
(fi', T") = (0, 0.4), Fig. 1 shows there are two degenerate 
global minima in $ at n = n = and n = n = 1; they 
give respectively {xpip) / (2N c ) = —1 and 1, and change 
from one to another through a chiral transformation or 
(n, n) (1 - n, 1 - ra). At (//, T') = (3, 0.4), Fig. 2 tells 
there are three degenerate global minima in $ at (n, n) = 
(0.1448, 0), (1, 0.8552), and (0.9519, 0.0481); they lead to 
(ipip)/(2N c ) = -0.8552, 0.8552, and respectively. At 
(fi',T') = (3.5,0.4), Fig. 3 indicates there is only one 
global minimum in $ at {n,n) = (0.9809,0.0191); this 
leads to (tpip) / (2N C ) = 0. 



III. PHASE STRUCTURE ON THE (fx, T) PLANE 

At the second order phase transition where the chiral 
condensate and the dynamical mass of quark vanish con- 
tinuously, there is only one global minimum in $, where 



d&/dn = d&/dn = 0. Using these conditions, we obtain 
an equation for the critical line 



= 1 



r 2 ) 



l + 3r 2 



T^ln- 



1 



1 



1 T 'a 
l+3r 2 



l+3r 2 



(15) 



which is depicted by the dotted line for r = 1 in Fig. 4. 
In the lower and left corner, the chiral condensate and 
dynamical mass of quark are non-zero. In the other side, 
they vanish identically. 

Below some finite Tg, the situation is different. There 
is a first order chiral phase transition line 



1 + 2^ 



(16) 



from some finite T% down to T 



0. this is illustrated by 
the solid line for r = 1 in Fig. 4. Chiral condensate and 
dynamical mass of quark jump from non-zero for // < ji' c 
to zero for [i' > [j! c . This is consistent with the data ob- 
tained by minimizing $, as shown in Figs. 1, 2 and 3. [i! c 
is larger than the dynamical mass of fermion at zero tem- 
perature. The reason, as explained in details in Ref. [12], 
is due to the fact that Wilson fermions break explicitly 
the chiral symmetry This seems counter-intuitive, since 
from a thermodynamical point of view, a transition is ex- 
pected when the chemical potential equals the dynamical 
fermion mass. It would be interesting to see whether the 
difference disappears in the continuum limit. 

The points when two lines described by Eq. (15) and 
Eq. (16) join are (n' 3 ,T^) and {/J,' 3 ,T"). (fj,' 3 ,T 3 ) is the tri- 
critical point, while the other one at high temperature is 
just a critical point on the second order phase transition 
line. Table I gives the location of the tricritical point for 
various r. For r = 1, we find {/J.' 3 ,T^ = (3,0.4498), i.e. 
the circle in Fig. 4. The phase structure for any r ^ is 
qualitatively the same. Note that the system along the 
line n' = ji' 3 + + experiences an inverted behavior: For 
T' € [0,Tg), the system is in the chiral-symmetric phase; 
While for T G (T 3 , T"), the system enters into the chiral- 
broken phase; For T" > T", the system is again in the 
chiral-symmctric phase. The value of T" is also given in 
Tab. I. Such a behavior exists even for naive fermions, 
though there is no tricritical point at finite T when r = 0. 

Figure 5 shows the results for the chiral condensate and 
quark number density (14) as a function // for T" = 1 
and r = 1, above the tricritical point; The results in- 
dicate there is a second order chiral phase transition at 
p! = 3.1770. Figure 6 shows those for V = 0.25 and 
r = 1, below the tricritical point; There is a first order 
chiral phase transition at \j! = 3. Figure 7 shows the 
chiral condensate for fi' = and r = 1; There is a second 
order chiral phase transition at T" = 2. These results are 
obtained by locating the minimum of $, and the critical 
point is consistent with the prediction of Eqs. (15) and 
(16). 
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IV. DISCUSSIONS 



In the preceding sections, we have investigated the 
QCD phase diagram in Hamiltonian lattice formulation 
with Wilson fermions. At the strong coupling, we find 
a tricritical point on the (n,T) plane, which has not 
been found in previous work in the Hamiltonian for- 
malism with Kogut-Susskind fermions or naive fermions. 
Our findings imply that on the (/z, T) phase, the phase 
structure of QCD with Wilson fermions (without species 
doubling) might be qualitatively different from naive or 
Kogut-Susskind fermions (with species doubling). Fur- 
ther detailed lattice study will be very important for un- 
derstanding the QCD phase diagram. 
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the system is in the chiral-broken phase. 
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FIG. 2. The same as Fig. 1, but for r = 1, fi' = 3.0, and 
T' — 0.4, where the system is at the criticality. 



r = = 3.5, T 1 = 0.4 
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FIG. 3. The same as Fig. 1, but for r = 1, fi' = 3.5, and 
T" = 0.4, where the system is in the chiral-symmetric phase. 
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FIG. 4. Phase diagram. The solid and dotted lines stand 
respectively for the first and second order transitions. The 
circle is the tricritical point. 
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FIG. 6. The same as Fig. 5, but for r = I and T' = 0.25. 
Eq. (16) predicts a first order chiral phase transition at 
He = 3. 
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TABLE I. Dependence of tricritical point on r. The point 
T" described in the text is also given. 
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FIG. 5. Chiral condensate and quark number density as a 
function of // at r = 1 and T' = 1. Eq. (15) predicts a second 
order chiral phase transition at fi'c = 3.1770. 
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FIG. 7. Chiral condensate as a function of T' at r = 1 
and \j! = 0- Eq. (15) predicts a second order chiral phase 
transition at T' c = 2. 



